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Abstract 

We reconsider the importance of statistical fluctuations for fission dynamics beyond the 
saddle in the light of recent evaluations of transport coefficients for average motion. The 
size of these fluctuations are estimated by means of the Kramers-Ingold solution for the 
inverted oscillator, which allows for an inclusion of quantum effects. 
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Kramers' transport equation jl] delivers the classic description of the dynamics of fission 
at finite excitation. It restricts to the high temperature regime in which collective motion is 
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treated within the frame of classical statistical mechanics. An extension to the quantum case 
has been described in [0] , and ]4| . This is possible within a locally harmonic approximation, 
in which global motion is described in terms of local propagators. The latter are constructed 
as special solutions of an appropriate transport equation in which anharmonic forces are 
linearized around some given point in phase space ||. For the present purpose this equation 
is needed only for the region between barrier and scission point . To simplify matters we like to 
apply the schematic model of ||, in which the potential is represented by an inverted oscillator 
(centered at q — 0) and where the coefficients for inertia M and friction 7 do not change along 
the fission path. The transport equation for the distribution d(q, P, t) in collective phase space 
may then be written as: 



d_P_ d_ d_ P_ d 2 d 2 

d~q~M + dP q+ dP^M + qp dqW + pp dPdP 



d(q,P,t) 



with C being the (negative) stiffness coefficient. There is little doubt that a correct treatment 
would require to go beyond the harmonic approximation underlying the form ([!]). Not only 
that the potential energy along the fission path will, in general, be more complicated than 
given by the simple quadratic dependence Cq 2 /2. Also the other transport coefficients will 
vary with q. However, recent numerical computations 0, || have revealed that the following 
ratios 

are quite stable. As we shall see soon, within the harmonic approximation it is them which 
parameterize the quantities we are mostly interested in, namely the kinetic energy and its 
fluctuation. 

As compared to the form of Kramers' equation, only two modifications occur in (|IJ), both 
referring to the diffusive terms. Firstly, there appears a cross term with the D qp being different 
from zero in the quantal regime. Secondly, the coefficient D pp will be given by the classic 
Einstein relation D pp = 7T only at large temperatures when quantum effects disappear. In 
their quantum version, these coefficients are defined by the following expressions 

^ X>K^SS (3) 
D qp = i^h coth (^) *» \C - Mu 2 } = C^ q - ±Z; p (4) 
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They are a consequence of the fluctuation dissipation theorem (FDT): Together with some 
simple symmetry relations, it allows one to calculate the "equilibrium fluctuations" E^j and 
E^ from the dissipative part Xqqiu) of the response function which represents average motion 
in q(t). For the linearized version used here, the latter must be related to the one of a damped 
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oscillator determined by M, 7 and C. If (0) is applied to a bound oscillator with C > 0, the 
expressions (|3], ^) warrant that for t — > 00 the dynamical fluctuations in g, P turn into those 
of equilibrium as determined by the FDT 0. Thus the contour C has to be chosen in the 
common way, namely along the real axis extending from —00 to 00. For unbound motion with 
C < 0, the case we want to study here, one has to apply suitable analytical continuations. 
This is possible in two ways: (i) One may evaluate the forms given in ([| |j) for C > and 
perform the continuations in the expressions one obtains after performing the integration (see 
[[| and (U). (ii) In an alternative method [|l(J one redefines the contour C. This allows for 
more general applications like in |TT[| . For more details we refer to |§ and ||. 

With respect to the evaluation of (|3], ^) we need to clarify a problem hidden in the integrals. 
In the quantum regime, some of them would diverge if we one were to take for the response 
function that of the damped oscillator. This problem is well known, and one possible solution 
is to apply the Drude regularization. This means to replace the Xqqi^) by 



In this way a "cut-off" frequency ud is introduced. We do not want to discuss the interesting 
questions of how its value can be fixed and from which physical quantities. For the compu- 
tations to be discussed below we chose ojd = 10zu, with the w given by @. Fortunately, the 
diffusion coefficients do not depend on ujd too much (see ||). Changing the latter by a factor 
of two, our final results would have to be modified by less than 30% which, as we shall see, 
will not influence much the conclusions we are going to draw below. 

Notice, please, that this regularization problem disappears in the classical limit. The 
latter is obtained if the h coth (hu / (2T)) is replaced by 2T/u>. With such a weighting factor 
all integrals in (||) and (f|) converge even for the case of a constant friction force. Apparently 
this classical limit is identical to the high temperature limit, for which one needs to have 
Tiu <C T. Looking back to the right hand sides of (Q) and (f|), it becomes evident why in this 
case the diffusion coefficients turn into those of Kramers' equation. For stable modes this is 
simply a consequence of the (classical) equipartition theorem. Since the diffusion coefficients 
then become independent of C the analytic continuation is trivial, and does not change the 
results for D pp and D qp when turning to the unstable situation. 

Whereas for stable modes the extension of Kramers' equation to the quantum regime is 
possible for all temperatures and all possible values of the transport coefficients, for unstable 
ones the description ceases to make sense at very low temperatures. First of all, below a 




(5) 



with a frequency dependent friction coefficient 




(6) 
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certain To it is not possible anymore to save the integral representation of the FDT. The 
contour C needs to cross the imaginary axis between the pole uj + = i \ uj + | of the unstable 
mode and the first Matsubara frequency oj\ 4 = 2nT/h (see Pl-lLl]). The T obtained in 
this way, namely T = hu + | /{2tt), is identical to the so called "cross over" temperature 
known from treatments of "dissipative tunneling" with functional integrals in the imaginary 
time domain. Here we are looking at real time propagation for which in a model case Ingold 
[ p"^| has been able to construct a phase space distribution corresponding to a constant flux 
across the barrier. Such a construction is possible above a critical temperature T c > T . As 
shown in || this distribution solves the transport equation ([!]) if only the diffusion coefficients 
are defined as described above. In this sense it may be considered the generalization to the 
quantum regime of the stationary solution found by Kramers. Actually, this T c turns out to 
be that temperature at which the CS^ becomes negative. It is of the order of 0.5 MeV or 
less, depends on the transport coefficients of average motion and decreases with increasing 
damping, see ||. On the nuclear scale such values of T can be considered small. (As a 
matter of fact, in such a regime the very concept of temperature itself becomes questionable.) 
Commonly, the fission experiments, which one may want to interpret with such a transport 
equation, involve higher excitations, for which our extension thus applies. Unfortunately, an 
experimental verification of the existence of quantum effects is still missing. Two possibilities 
have been suggested so far, the (dwell) time r from saddle to scission || (see below), and 
the decay rate 0. For those, however, quantum effects would show up only in a very narrow 
range between T c and values of about T ~ 1 • • • 1.5 MeV. 

This situation may change if one looks at quantities which involve the momentum distri- 
bution. Such a feature is known from studies of the dynamics of stable modes. There the 
equilibrium fluctuations in the coordinate get squeezed when friction increases. In this way 
their values get closer to the classical limit. The opposite holds true for the momentum. A 
nice demonstration of this effect can be found in |H| where path integrals are applied to a 
solvable model. In |§ this problem has been taken up for the nuclear context within the 
locally harmonic approximation, for which one is not restricted to describe the "heat bath" 
of the nucleonic degrees of freedom by a set of coupled oscillators. 

In the present letter we specifically address the dynamics across the fission barrier. Dif- 
ferent to and m we want to exploit the Kramers-Ingold solution di(q, P) to evaluate the 
kinetic energy and its variance at scission. The calculation can be done in complete analogy 
to the case discussed in || for Kramers' equation. The average kinetic energy at scission may 
be defined by 

r°° PdP P 2 
Bm„=/ — d^F) — (7) 
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and its variance by 



a 



kin 



°° PdP . , ^ ( P 2 \ ,„ , 2 



with the g sc being the coordinate of the scission point. Notice please, that the sampling is 
done for the following normalization of the distribution 

00 PdP , . . 

di(g,P) = l (9) 



00 



with j being the constant flux across the barrier. Like in the classical case (see eqs.(13) and 
(14) of ||) the integrals can be carried out analytically. Introducing the abbreviation 



the result writes as 



and 

a 2 



C = jl + rP- V = (10) 

£kin = \ [Mw\l c - Dgpj C 2 + 5* (1 + (11) 




(12) 



For the difference AV = \Mw 2 q 2 c of the potential energy between saddle and scission the 
value AV = 20 MeV will be adopted in the following. Notice please, that the first term in 
(|TID stands for the kinetic energy = Mzu 2 q 2 c ( 2 /2 one would obtain in a mere trajectory 
calculation. This is easily verified from the solution 

P 1 

qsc = g(T) ~ 2M^+^ eXp(U7( V 1 + ^ ~ ^ )r) (13) 

of Newton's equation for q(t = 0) = 0, only assuming the saddle to scission time r to be 
sufncently large to have zjy/l + i] 2 t ^> 1. In (|TTD all the other terms represent the effects of 
the fluctuating force. 

Before turning to discuss numerical evaluations of the expressions just presented it may 
be worth while to comment on their physical relevance. Firstly, we like to stress that by using 
the stationary solution dj(q, P) the results become insensitive to initial conditions, like those 
on top of the barrier one would need to invoke in a time dependent picture. Please notice 
that the only uncertainty left in di(q,P) is the multiplicative factor hidden in the current j 
which drops out when calculating the -Ekin and < Jy m according to (0) and @, respectively. For 
the situation to which this stationary solution dj(q, P) of the inverted oscillator commonly is 
applied to (c.f. 0), this j can be said to stand for the decay rate out of the potential miminum. 
Eventually, one would then like to have this minimum to be well pronounced, in the sense 
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of having the barrier height be large compared to T, more precisely to that temperature 
the system has on its way towards the saddle. However, as demonstrated in [14 and |L5| 
the di(q, P) may as well be understood to result from integrating a t-dependent distribution 
d(q, P, t) over time t. It is not difficult to understand that such a time integrated function 
solves an equation like ([!]) with the left hand side put equal to zero. Moreover, as shown in 
[ [0| , already for quite small damping rates 77 (called 7 in [15]) this time-integrated distribution 
shows a relaxational behavior to Kramers solution, if considered in its dependence on q. In 
this spirit, formulas (0) to fll~2| ) may be understood in the following way. Provided that the 
scission point q sc does not lie too closely to the position of the barrier top, we may assume 
the di(q,P) to adequately portray the momentum distribution of an actual fission process. 
Formulas (|ll|) and ([12]) then measure the kinetic energy and its fluctuation for the distribution 
one obtains after summing up all events leading to fission. 

In Fig.l we show the kinetic energy as function of 77 for different temperatures as calculated 
from (|ll|)- Fully drawn lines correspond to the quantal case and dashed ones to the high 
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Figure 1: Average kinetic energy as function of 77 for different temperatures T. Fully drawn 
curves include quantum effects, dashed ones refer to Kramers' solution. Shaded dots give 
results based on the values for w and 77 (at given T) as obtained in microscopic computations 
of [7] (dark shaded) and [8] (light shaded). (In this and the following figures we use units of 
h = 1). 
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temperature limit, thus corresponding to Kramers' equation. As the figure demonstrates, 
quantum effects increase with damping and they amplify the average kinetic energy. As a 
matter of fact, the latter is seen to attain quite large values in any case, if compared to the 
value of AV. Obviously, this is an effect of the fluctuating force, as friction alone acts to 
diminish the velocity and thus the kinetic energy. This feature is demonstrated explicitly in 
Fig.2. There the ratio E^/E^ is shown (with the E^ } introduced below fli~2l)). In both 
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Figure 2: Ratio of the kinetic energies, E^ and E^ n , calculated, respectively, without and 
with fluctuating force. The dots are specified as in Fig.l. 

figures shaded circles indicate results obtained by using the transport coefficients for average 
motion as found in microscopic computations of [0] (dark shaded) and (light shaded). 
Fig.2 demonstrates clearly that trajectory calculations may grossly underestimate the size of 

E-kin. 

The very fact of the big influence of the fluctuating force hints at the importance of 
statistical fluctuations of the kinetic energy itself, which may be calculated according to ([12]) . 
In Fig. 3 we show the square root of the variance a divided by the kinetic energy -Ekin, both 
calculated at the scission point. For large damping and large temperatures this ratio comes 
close to the limiting value y2/3. When the effective damping rate r\ is somewhat larger than 
1 this value of y 2/3 is reached for practically all T. Incidentally, we may note that this 
ratio cr/E'km can be expected less sensitive to the Drude frequency ojd than the individual 
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Figure 3: Variance a over average kinetic energy Ey m as function of T/zu, for different values 
of r\. 

quantities, simply because both of the quantities, a and £?kin, change with u D alike. 

Next we like to work out more explicitly the size of quantum effects. In Fig.4 we plot for 
both the kinetic energy and its fluctuations, the ratio of the values in the quantal case to the 
corresponding ones in the high-T limit (Kramers' case). As seen from the figure, these ratios 
may take on quite large values and they increase with increasing damping. Thus Fig.4 agrees 
with the observation made in Fig.l and confirms our conjecture raised earlier: Whenever the 
collective momentum is involved quantum effects get larger with increasing damping. Like 
before we have again indicated by dots the range one would expect for these values on the 
basis of the microscopic computations of [0 and |§ . Notice, please, that in these computations 
the coefficients r\ and w have been evaluated as function of T. For the curves shown in all the 
figures these coefficients have been varied as free parameters. This will facilitate comparison 
with other theoretical models and with results deduced from experiments. Indeed, it may 



be said that there is experimental evidence (see e.g. |16| and ||17|| ) for much larger damping 
coefficients than found in [0] and ||. The authors of [[H]] need values of as much as rj = 10 to 
cope with findings in experiments where fission of heavy nuclei is observed accompanied by 
GDR 7 rays, at temperatures not larger than about 2MeV. Applying macroscopic pictures 
to evaluate the transport coefficients, like the wall formula for friction, irrotational flow for 
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Figure 4: Ratios of kinetic energy and its variance to their values in the high temperature 
limit, both as functions of T/w, again for different 77. 



inertia and the liquid drop model for the stiffness one would find similar values [H| (for 224 Th 
at T = 2 MeV r] becomes ~ 7.5). If one looks at lighter systems where large angular momenta 
are needed to find fission events with some finite chance , the barrier becomes quite small 
and broad, thus leading to smaller values of w and hence to larger values of rj As seen 
from our figures, this hints not only at the importance of the statistical fluctuations as such, 
but also at the necessity of calculating them with the quantal diffusion coefficients. Most 
likely this may again modify the interpretation of experimental results. In any case, it is 
probably fair to say that still some work is to be done before more conclusive statements 
can be made about the size of the transport coefficients. In this context one may mention 
that there are indications from other experiments that rj might be smaller, indeed, than the 



macroscopic picture requires 20 



So far in this paper we have been looking at cases which explicitly involve the momentum 
distribution in one way or other. For the sake of completeness we should like to take up once 
more the question of the influence of statistical fluctuations on the time r it takes for the 
system to move from saddle to scission. In || the following formula had been derived, based 
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on Kramers' equation: 



VI + V 2 - V + 7]' 2 - 7]) 

with lZ(x) being the Rosser function 

rx poo 

1Z{x) = I exp(y 2 )dy / exp(— z 2 )dz 
Jo Jy 



(14) 



(15) 



-\Cq 
-*qq 



In the quantum case the argument of the 1Z in (|14]) would get the additional factor yT/CS 
[H], but as can be seen from Fig.2 of this reference this modification may safely be neglected 
for T/hw > 0.5. For that regime Fig. 5 shows the wt as function of AV/T = Mw 2 q 2 c /2T. 
As the hw is of the order of 1 MeV, the Tq is seen to lie in the range of 1 — 2 x 10~ 21 s. 
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Figure 5: The mean saddle-to-scission time tq (scaled by to) as function of the difference 
AV of the potential energy (measured in units of T) . 



It is interesting to compare the result (|TJ]) with the one one would get if the system would 
move along a trajectory starting on top of the barrier with the momentum Pq (i.e. including 
disipation but discarding the fluctuating force). As seen from (|13|) one gets a formula like (|i~4l) 
but with the 21Z(^J Mw 2 q 2 sc /2T) replaced by hx(2Mwq sc ^/l + r] 2 /P ). To make the analogy 
to ( |i~4|) even closer one may estimate the initial momentum by associating P 2 /2M to an 



average, thermal kinetic energy on top of the barrier. One might be tempted to use for the 
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latter E^ n = T/2, the value given by the equipartition theorem. A numerical evaluation 
shows that in this case one would overestimate the r by about 10% (for 77 = 1) to 50 % 
(for i] = 10). However, as we may learn from (|TT|) , for the inverted oscillator the stationary 
solution suggests a larger value of the average kinetic energy. Putting there the q sc = one 
obtains (in the high-T limit) E^ n = T(l + r](), which for zero damping gives twice the value 
of the equipartition theorem. This modification would improve the results slightly. 

Finally, we wish to add some remarks on the variation with q. It is evident from the 
discussion given above that under certain circumstances the Kramers-Ingold solution may 
be applied for regions before scission. It may thus be used for performing averages over 
quantities which simply depend on the collective variable. Prime examples are the evaporation 
probabilities for light particles |2l| and 7-rays (see [171). Adopting the same normalization 
as in (|7|, H) one gets as weighting factor the current j which according to @ is constant. 
This implies that for the motion from saddle to scission an average in the coordinate q simply 
reduces to the algebraic one. 
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